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ABSTRACT 

We study the evolution of baryonic gas before the reionization in the 
lognormal (LN) model of cosmic clustering. We show that the thermal history 
of the universe around the reionization can roughly be divided into three 
epochs: 1) cold dark age z > Zr, in which baryon gas is neutral, and opaque to 
Lya photons; 2) hot dark age Zr > z > Zgp, in which a predominant part of 
baryon gas is ionized and hot, but it is still opaque to Lya photons; 3) bright 
age z < Zgp, in which the universe is ionized highly enough to be transparent 
to Lya photons. In the flat cold dark matter cosmological models given by 
WMAP and COBE, the difference of the two redshifts found to be 

as large as ~ 10 with Zr ^ 17 and Zgp ~ 7. This reionization history naturally 
yields a high optical depth to the CMB Te ~ 0.12 - 0.19 observed by the TE 
polarization of the WMAP, and a low redshift Zgp of the appearance of the Lya 
Gunn- Peterson trough 2;pp ~ 6 — 8 in QSO's absorption spectra. The reason 
why the universe stays long in an ionized, yet Lya opaque, stage is because the 
first photo-ionization heats the intergalactic gas effectively and has balanced the 
gravitational clustering a long period of time. Therefore, the result of a high 
Te and low Zgp is a common feature of all the models considered. Besides the 
cosmological parameters, the only free parameter we used in the calculation is 
Nion, the mean ionization photons produced by each baryon in collapsed objects. 
We take it to be 40 - 80 in the calculation. 

Subject headings: cosmology: theory - large-scale structure of the universe 
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1. Introduction 

The large scale structures of the universe evolved from initial mass density 
perturbations seeded by the quantum fluctuations of inflation. The power spectrum of the 
initial perturbation is of power law, P{k) oc A;", with index n ~ 1, and so the mass density 
perturbation is stronger on smaller scales. Accordingly, objects on smaller scales generally 
formed at earlier time. Small objects like primordial black holes could form at the early 
radiation-dominated era. Objects on scales of stars could form in the dark age, i.e. between 
the recombination and the reionization. These stars play the leading role of the history of 
the reionization. 

The study of structure formation in the cosmic dark age is challenged recently by two 
observations. 1. The complete Lya Gunn- Peterson trough in the absorption spectra of high 
redshift quasars shows that the reionization occurred probably at the redshift 2:^ — 6 — 8 
(Fan et al. 2002, and references therein). 2. The TE polarization of the Wilkinson 
Microwave Anisotropy Probe (WMAP) yields a model-independent value of the electron 
scattering optical depth to the cosmic microwave background (CMB), Te = 0.17 ± 0.04, 
which requires the reionization redshift to be 2; ~ 11 — 30 if the reionization happen rapidly 
(Kogut et al. 2003). Therefore, a problem naturally arisen is how to reconcile the two 
results, or why the reionization redshift from the Lya spectrum is significantly lower than 
that given by the optical depth to the CMB. It has been realized that the high Te and 
low Zr probably indicates that the reionization is not abrupt, but lasts a long period from 
z = 30 - 20 to 6 - 8 (Cen 2003; Haiman & Holder 2003; Whyithe & Loeb 2003; Holder et 
al. 2003; Hui & Haiman 2003; Onken & Miralda-Escude 2003). The further questions we 
might ask are how to describe the long period of the reionization and which is the physical 
mechanism leading to such a long period of the reionization? This paper is to address these 
problems using a semi- analytical approach in the lognormal (LN) model. 
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The LN model assumes that the mass field of intergalactic medium (IGM) is given 
by an exponential mapping from the linear field (Bi 1993; Bi & Davidsen 1997). The 
probability distribution function (PDF) of the IGM field is then lognormal. The first 
argument in supporting the lognormal PDF is probably from the isothermal model of the 
IGM in the gravitational potential as 

(1) 

where p(x) is the IGM mass density and (f)dm{^) the gravitational potential. T and n are, 
respectively, the temperature and molecular weight of the gas, and 7 is the polytropic 
index. Eq. (1) is the well-known IGM model of intracluster gas (Sarazin & Bahcall 1977). 
Zeldovich et al. (1990) pointed that when the gravity potential 0d,n(x) is a random field 
with a Gaussian PDF, the field p(x) of eq.(l) is lognormal. It has been shown that the 
lognormal solution is a reasonable approximation of the dynamical equations of the IGM 
(e.g. Jones 1999; Matarrese & Mohayaee 2002). Therefore, the LN is dynamically legitimate 
to model the weakly nonhnear properties (Bi & Davidsen 1997), and highly non-Gaussian 
features (Feng & Fang 2000) of the IGM. Recently it is also used to study the collapsing of 
hydrogen clouds of the IGM (Bi et al. 2003, here after BFFJ). 

As has been emphasized in BFFJ, the formation of the first generation of collapsed 
hydrogen clouds is given by rare events, which is sensitive to the tail of the mass field PDF 
on the high perturbation side. It is noted that the LN PDF is of long tailed and predicts 
a higher probabihty of rare events of high perturbation than PDFs without a long tail. 
Upon this consideration, the BFFJ showed that the reionization lasts a period from redshift 
15 to about 10. However, the BFFJ did not considered the reaction of the evolved IGM 
on the reionization. During the reionization, the IGM underwent a dramatic evolution 
from neutral to ionized hydrogen, from cold to hot gas, and from low to high entropy 
states. Consequently, the clustering behavior of the IGM should also undergo a remarkable 
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evolution in the epoch around the reionization. On the other hand, the clustering of the 
underlying dark matter is not affected by the reionization. Therefore, we must analyze the 
impact of the IGM evolution upon the reionization history. 

The LN model directly provides the spatial distributions of the IGM. It is effective to 
reveal the reaction of the evolution of the IGM on the reionization. We will show that a 
long period of the recombination is a natural outcome of the LN model if considering the 
evolution of the IGM at the reionization. A long period of the reionization history generally 
yield a high optical depth Tg to the CMB and a low reionization redshift to the Lya forests. 
We will calculate these quantities in the cold dark matter cosmological models. 

The paper is outlined as follow. In §2 we will describe the LN model of the reionization. 
The history of hydrogen reionization will be discussed in §3. Finally §4 gives the discussion 
and conclusion. 

2. The Lognormal Model 

2.1. Mass Field of Baryonic Gas 

In the LN model, the mass field of the IGM is given by an exponential mapping of the 
linear density field of the underlying dark matter as 



where (Tq is the variance of the linear Gaussian field of the dark matter filtered on scales of 
the Jeans length of the IGM; Sq{x.) is the linear density contrast of dark matter ^^^(x) 
smoothed by a window on the Jeans scale, i.e. its Fourier counterpart is (Bi & Davidsen 



p(x) = pexp[5o(x) - ao/2]. 



(2) 
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where Xf, is l/27r of the comoving Jeans length of the IGM. 

Obviously, when the fluctuation (5o(x) is small, i.e. (5o(x) <^ 1, eq.(2) reduces to 
approximately 5(x) ~ 5o(x), where 5(x) = [p(x) — is the density contrast of the IGM. 
Moreover, when the fluctuation 5o(x) is large, combining cqs.(2) and (3) yields eq.(l) (Bi 
& Davidsen 1997). Since the initial density perturbation is stronger on smaller scales, the 
perturbations of the dark matter gravitational potential 0cim(x) are larger on smaller scales. 
Thus, even at early universe the IGM field on small scales may already be in the state 
described by eq.(l). Thus, it would be not appropriate to assume that the IGM field is 
always Gaussian on all scales. The LN model provides a uniform description of the IGM 
distribution in the linear, weakly nonlinear and lognormal regime at various scales and 
redshifts. 

Since the random field 5o(x) is Gaussian, the PDF of the field p(x) is lognormal as 



p{p/p) = z ^= cxp 

(p/p)(To\/27r 



1 (\n{p/p) + al/2 



2 



, p>0. (4) 



2 V cfo 

which gives the probability density of the event such that the mass density within an area 
on scales of the Jeans length is found to be equal to ~ p/p. The tail of the PDF eq.(4) is 
~ exp[— (l/2)(lnp/p)^/(7Q + Inp/p]. Therefore, the probability of p/p for a given ctq would 
be equal to the probability of {p/pY for 2(7o. Namely, the change of variance ctq by a factor 
2 leads to a change of the possible high density events from p/p to {pjpf. 

Obviously, the PDF is normalized as p{x)dx — 1, and the variance of p/p is given 
by (Vanmarcke 1983) 

a = [e'^l - (5) 



The n*'* moment of p is 



{p/pY = exp 



(n^-n)^ 
^ ' 2 



(6) 



Equation (6) shows [(p/p"]^/" > [{p/ pYY^"^- That is, the ratio between the high- and 
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low-order moments is divergent in the limit of n — > oo. A field with such a divergent 
moment ration, by definition, is intermittent (Zel'dovich, Ruzmaikin, & Sokoloff 1990). 
Actually, LN is typically a mathematical model of intermittent field. This property is 
in good consistent with the detected intermittency of the transmitted flux of QSO Lya 
absorption (Jamkhedkar, Zhan & Fang 2000; Feng, Pando & Fang 2003, Jamkhedkar et al. 
2003). 



2.2. Jeans Length of Baryonic Matter 

Equation (2) shows that, in the LN model, the clustering features of the baryonic gas is 
characterized by the variance uo of the linear density field on the scale of the Jeans length. 
Primordial baryons, created at the time of nucleosynthesis, recombines with electrons to 
become neutral gas a.t z ^ 1000. Thereafter, the gas cools down adiabatically with the 
expansion of the universe. In a homogeneous universe with the mean mass density pm, the 
Jean length of gaseous baryonic matter or the IGM is defined by A;, = Vs{n/Gpmy^'^, where 
Vs is the sound speed of gas. The corresponding Jeans mass is mj = (47r/3)Afpm- In a 
comoving scale of = Xb/^ir, we have (Bi & Davidsen 1997) 

nl/2 



2-fkBT 



(7) 
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The polytropic index is assumed to be 7 = 5/3, and hydrogen temperature thus follows 
T oc p^/^, where p is the mean mass density of baryonic matter. At 2; = 15, the hydrogen 
temperature is ~ 4 K (Medvigy & Loeb 2001), and the Jeans length Xb ~ 1.2 kpc. The 
redshift evolution of the comoving Jeans length Xb is approximately given by (1 + zY^'^. 

During the reionization, the baryonic gas will be heated by the UV ionizing background 
from few degrees of temperature to ~ 1.3 x W^K; and correspondingly, there is about 2 
magnitudes increasing in Xb or 6 magnitudes in the Jeans mass (e.g. Ostriker & Gnedin 



-8- 



1996). After the reionization, the IGM temperature is maintained at about ~ lO^X by the 
UV background photons, and therefore, Xb will gradually increase with the decrease of z 
due to the factor {1 + z) in eq.(7). 

For the variance in cq.(2), we consider three fiducial cosmological models. Two are 
inferred from the WMAP, and one from the COBE. The cosmological parameters of these 
models are hsted in Table 1. We use the hnear power spectrum of the dark matter given 
by Eisenstein & Hu (1999). The linear variance is approximately growing with redshift as 
(To{z) oc 1/(1 + 2;). For cosmological models of Q = 1 and ^ 0, the linear growth of 
density perturbations is not exactly following the factor 1/(1 + 2;) (Lahav et al. 1991). We 
take account of this effect by replacing the factor 1/(1 + 2;) by l/g{Qrn: ^^a)(1 + z), where 
the g factor is 



This correction is approximately reasonable at 2; > 1 (Avehno & de Carvalho 1999, Carroll 
ct at. 1992). For the models of WMAPl, WMAP2 and COBE, the g factor is also shown 
in Table 1. 

Because the powers of the linear perturbations k^P{k) in all of these models are smaller 
on larger scales, the two order magnitude jump of xi, at the reionization generally leads to a 
drop-off of the variance (Jq. A typical curve of (Jq for the model WMAP2 is displayed in Fig. 
1, in which the reionization is assumed to happen instantly around redshift Zr = 17. The 
dotted line in Fig. 1 is for ctq in case of no reionization occurring. The zigzag 2;-dependence 
around the reionization is a common feature of the (Tq curve. For instance, the Fig. 2 of 
the BFFJ shows also a zigzag at 2; = 7. In the BFFJ, it is assumed that the reionization 
takes place at redshift Zr = 7. In this paper, we will not use a priori zigzag redshift Zr, 
but determine it by a self-consistent calculation. It should be pointed out that the zigzag 
2;-dependence of (Jq does not affect the evolution of dark matter, of which the clustering is 
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Table 1: Optical depth in three cosmological models 

model flm f^A fife Us (7s h g 

WMAPl 0.27 0.73 0.044 0.97 0.84 0.71 0.758 

WMAP2 0.29 0.71 0.047 0.99 0.9 0.72 0.772 

COBE 0.3 0.7 0.045 1 0.85 0.7 0.797 

basically regardless the reionization. 

2.3. Cumulative Mass Fraction and Clumping Factor 

For a field with the PDF of eq.(4), the cumulative mass fraction M(> p/p), i.e. the 
fraction of mass in regions having mass density larger than a given p/p, is 

POO I 

M[> {p/p)] = / xp{x)dx = -erfc 
Jp/p 2 

which implies that M[> {p/p)] is a strongly varying function of the variance (Jq, but weak 
for the variable p/p. Fig. 2 demonstrates a cumulative mass fraction for p/p = 2, 6 and 
10, in which ctq is taken from Fig. 1, and using the dotted line for z < Zr, i.e. without 
considering the drop-off of (Jq at Zr- Thus, Fig. 2 can be only used at redshift z > Zr- 

The decrease of M{> p/p) with p/p is much slower than a Gaussian PDF, or PDFs with 
a Gauss-like tail. For a Gaussian tail, we have M{> 10)/M{> 6) ~ exp[(-102 + 6'^)/2a^], 
or M{> 10)/M{> 6) < exp(— 30), when (Tq < 1. However, Fig. 2 shows that even when z is 
as high as 20, M{> 10) is less than M{> 6) only by a factor less than 10. This is obviously 
due to the long tail behavior of the lognormal PDF, which gives a larger probability of high 
density events at high redshifts when the variance ctq is still small. 

The variable p/p in eq.(9) is the density within an area on scales of the Jeans length, 
and therefore, the area with p/p > 1 has mass larger than the Jeans mass. However, these 
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regions are not always in collapsing. This is because that the IGM clouds have no time 
to collapse even when the cloud mass is larger than the Jeans mass. BFFJ found that at 
redshift 2; ~ 7, only clouds with density p/p>Q are in the turning phase. The threshold 
p/p ~ 6 is much larger than the Jeans mass. In other words, the clustering on the Jeans 
length scales is significantly suppressed. The suppression of collapsing on the Jeans length 
scales has also been noted in studying the so-called filtering scale (Gnedin & Hui 1998). 

The threshold for the clouds going to collapsing at redshift 2; > 7 is a httle lower than 
6. Therefore, the mass fraction of collapsed clouds can be estimated by the cumulative 
mass fraction M{> 6). One can estimated the velocity dispersion of the collapsed 
clouds by considering that in the LN model the collapsed clouds in potential well is 
probably isothermal [eq.(l)]. In this case, the velocity dispersion of the collapsed clouds is 
approximately (e.g. Jing & Fang 1994) 



Thus, for h = 0.7, Qm = 0.3, the velocity dispersion of the collapsed clouds of p/p > 6 is 
from > 3.8 km s~^ at ^ = 10 to > 7.7 km s~^ at ^ = 20. If using p/p > 10, we have > 5.0 
km at 2; = 10 to > 9.8 km s~^ at 2; = 20. Therefore, all these clouds have mass larger 
than 10^ Mq, the H2 coohng can take place. 

The so-called clumping factor is the second moment of the IGM density field defined by 



Fig. 3 plots the clumping factor C{z) in the three models with ao without considering the 
drop-off at Zr- It can only be used for redshift z > Zr- The clumping factor is significantly 
z-dependent. For all three models, C{z) are large only at redshift 2; < 10, and rapidly 
decreases to < 3 at redshifts z > 15. At a first glance, it seems to be a contradict between 
a higher probability of clustering and a lower clumping factor. Actually, it is nothing but 



a,c^l3HoXhil + zy/^nl/\p/py/^ 
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a typical feature of intermitted field, which consists of high density spikes and low density 
regions between the spikes, and therefore, dense objects can be effectively formed in the 
spikes, while the mean clumping factor over the whole space is relatively low. 

3. History of the Reionization 

3.1. Clustering of Baryonic Gas around Reionization 

From Fig. 1, one can see that the 2;-dependence of ctq before and after the reionization 
Zr are similar. Before Zr, (Tq is increased from 0.4 to 0.8, and after Zr, ctq is again to increase 
from about 0.4 to higher values. A fall of uo is by a factor of ~ 2 at 2;^. Of course, the 
change of ctq at the drop is dependent on Zr- In this section, we only qualitatively describe 
the effect of the (Jq drop on the clustering of baryonic gas. Quantitative calculation will be 
given in next two subsections. 

As emphasized in §2.1, the cumulative mass fraction is a strong function of ctq. The 
drop of (To by a factor of two leads to a substantial change of the cumulative mass fraction. 
One can then divide the IGM clustering into two phases, before and after redshift Zr- Before 
the reionization, structure formation was weak when the variance ~ 0.4, and became more 
active when the variance reached ~ 0.8. Just after (Jq dropping, the structure formation 
tended to inactive again as the variance reduced to ~ 0.4. Soon after the variance growing 
up to ~ 0.8, the structure formation was built up once more. Therefore, there are two 
phases, in which the rate of structure formations are high. The early phase is just before 
Zr, and the late phase is at the second time of (Jq > 0.8. 

The physical reason for the two-phase scenario of structure formation is obvious. 
Baryonic gas falls into dark matter gravity wells and the first generation of stars forms 
before Zr- Because the UV ionization photons background is gradually established from the 
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first generation stars, baryonic gas is heated to temperature 10 K, and the Jeans length 
becomes bigger. The irregularities or clustering within the big Jeans scale would be smeared 
out by thermal motion. Since the density fluctuations on the big Jeans scales would still 
be in the linear regime, the process of the structure formation is significantly suppressed 
even halted. At the end of this smearing out, some ionized hydrogen will move away from 
the bottom of the gravity wells of dark matter halos. Eq.(l) also shows that when the 
temperature T increases by a factor of 10^, many nonlinear structures will disappear. The 
clustering will recover till the variance ctq growing again. In the second phase, the objects 
form from ionized IGM, and their masses range in the order of 10^~^° M©. 

In the two-phase scenario, the structure formation will be slowed down just after the 
reionization. Consequently, there is no enough UV photons produced by stars to maintain 
a fully ionized IGM, and a significant fraction of HI remains in the IGM. In this case, the 
IGM is opaque for Lya photons. The IGM becomes transparent for Lyct photons till the 
redshift at which the rate of star formation grows again to be able to produce enough 
UV photons for maintaining a highly ionized IGM. Therefore, we can define two redshifts: 
(1) Zr, the redshift of reionization, which is the first time when the IGM was fully (or 
almost fully) ionized and heated up to temperature ~ 10^ K; (2) Zgp, the redshift of the 
Gunn-Peterson transparency, which is the starting time for IGM entering into the state 
such that the volume- averaged fraction of HI is less than 10^"^, and the Gunn-Peterson 
absorption ceased. In a word, the clustering of the IGM is suppressed during the period 
from Zr to Zgp, and thus leads to the reionization epoch lasting a long period from Zr to Zgp. 

The two-phase scenario is based on the z-ao relation indicated in Fig. 1. Obviously, 
this z-dependence of ctq is approximate, as we assumed that the temperature of the IGM 

suddenly increased by a factor of 10"'. More precisely, the sharp decline of ctq at Zr should 
be replaced by a soften transition from the top branch to the bottom branch in the z — ao 
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plane. Nevertheless, the mechanism for a long period reionization outlined above still works, 
as it essentially depends only on the decrease of ctq with the heating and ionizing IGM. 



3.2. The Redshift of Reionization 

Let us first calculate the reionization redshifts Zr- The reionization of hydrogen clouds 
is usually characterized by the filling factor of ionized hydrogen Qhii{z)- The reionization 
redshift Zr can approximately determined by requiring Qniiizr) > 1- The equation of 
Qhii{z) is given by (e.g. Barkana & Loeb 2001) 

dQmi NjQ^it) dM{> p/p) 1 

IT = It -sO^BCHiiitMHii, (12) 

where the factor 0.76 is from the fraction of hydrogen in the cosmic baryons, 
fiff — 1.88 X 10~^(n6h^/0.022) cm~^ is the present value of the mean number density of 
hydrogen, Nionit) is the mean ionization photons produced by each baryon in collapsed 
objects and might be t-dependent. The first term on the r.h.s. of eq.(12) accounts for the 
ionizing source. The ionization photons produced by primordial black holes are negligible 
because of the spectral index n < 1 on small scales (He & Fang, 2002). The second term 
on the r.h.s. of eq.(12) is from the recombination, in which a is the cosmic factor, is 
the case B recombination coefficient ~ 2.6 x 10~^^cm^s~^ for hydrogen at T = 10^ K, and 
Cniiit) is the clumping factor of HI I. 



Solving eq.(12) yields 



. . I^ NUt)dMi>p/p) 
^^^^^^^ - io ^?76 dF' ' ' ^^^^ 



where 

-0 /■* CHIljt") 

In the solution (13), an uncertain parameter is A^ionj which depends on the atomic processes 
in star formation. As usual, one can express this number by Nion — x fstar x fesc, where 
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fstar - the fraction of collapsed baryon that formed star; - the ionization photon number 
produced per baryon in star; and fesc - the fraction of photon that escaped from galaxy. 
Generally, the first generation stars are massive (Abel et al. 2002, Bromm et al. 2002), and 
their lifetime are much less than the Hubble time at high z that we are interested. Thus, it 
would be reasonable to approximate the photons being produced at the same redshift 
as collapsed clouds. 

The number is quite different for stars with different mass and metallicity. With 
the IMF of Scalo (1998), one has ?a 4000 for stars with metallicity ~ O.OSZ©. For the 
Salpeter's IMF, is ~ 6000 in average for stars with masses of 1 to 100 M0. For stars with 
zero metallicity, the number could be as high as 10^ (Bromm et al. 2001, Venkatesan 
et al. 2003). Observational estimations of fstar and f^sc generally yield f star fesc ~ 10~^. To 
sketch the reionization history driven by the mechanism discussed in §3.1, we use a simplest 
model of A^ion, i-e. let Ajon to be constant (^^-independent), and takes the value of 40, 60, 
and 80 respectively, which covers the possible range of Ajon- 

Figure 4 plots the fiUing factor Qhii vs. redshift z in the three cosmological models 
with Nion = 40, in which we take Chii{z) = C{z). The curve of the model WMAP2 is 
similar to the COBE, while the WMAPl shows a lower filling factor. This is because both 
the spectral index Ug and the normalization of power spectrum (t% of the WMAPl are 
smaller than other two models, and accordingly, the WMAPl gives smaller variance ctq 
on the Jeans length of the IGM. From Fig. 4, the redshift of reionization is in the range 

Zr> 16. 

We also estimate Zr by taking A^jon=40, 60, and 80, and p/p = Q, 8, 10. The results are 
tabulated in Table 2. In all the considered parameter-range, the redshift of reionization 
Zr is > 10. This result is consistent with the BFFJ, which also found Zr > 10 but using 
different method. The BFFJ's conclusion is based on the evolution of the mass density and 
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Table 2: Reionization redshift (zr) 



p/p 


6 


6 


6 


8 


10 


TV- 

'ion 


40 


60 


80 


40 


40 


WMAPl 


15.8 


17.0 


17.4 


14.0 


12.9 


WMAP2 


18.0 


19.5 


20.4 


16.0 


14.9 


COBE 


17.8 


19.0 


20.0 


15.9 


14.6 



velocity profiles of hydrogen clouds. Clearly, both approaches arrives at the same result 
that the reionization could occur much earlier than the redshift of an appearance of the 
complete Gunn-Peterson trough of Lya absorption. 

The assumption of Chii{z) — C{z) is correct only if the density distribution of HI I 
follows the distribution of total IGM. Since HII may favor to locate nearby collapsed 
clouds, we should use a biased relation Chii{z) = bC{z). The bias parameter b can be 
estimated by 



where we used eq.(6). Since 0.4 < ctq < 0.8 (Fig. 1), we have 2 < 6 < 5. These values of 
b is consistent with the estimation by the property of collapsed clouds (Mo & White 1996, 
Fang & Jing 1998). Because the velocity dispersion of collapsed clouds is small [Eq.(lO)], 
the parameter b should be in the above-mentioned range. Figure 5 plots the fiUing factor 
Qhii vs. redshift z for the model WMAP2, pjp — Q and Nion — 40, but taking the 
biased clumping factor of HII with parameter b = 1.5, 2 and 4. The results show that 
the uncertainty given by the bias factor is about the same as the uncertainty given by 
cosmological parameters (Fig. 4). Therefore this uncertainty will not affect the conclusions 
shown below. 
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3.3. The Redshift of the Gunn-Peterson Transparency 

Although Qhii is ~ 1 at 2;^, it does not mean that Qhii — 1 for all z < Zr, because Uq 
fell at z = Zr- For a self-consistent calculation of Qhii in the epoch of z < Zr, we should use 
the lower branch of ctq as shown in Fig. 1. 

The filling factor Qhii at z < Zr can also be calculated using equation (12). Figure 6 
plots the results of Qhii vs. redshift in the model WMAP2 with the parameters {p/p, Nion) 
of (6, 40), (6, 80), and (8, 40). Other models yield similar curves of Qhii- A common 
feature is that the reionization Qhii drops to a level less than 1 after Zr, and then keeps 
~ 0.5 in the epoch of redshift Zr > z > Zgp. It is due to the decrease of the star formation 
rate when ctq drops back to ~ 0.4. When z — > Zgp, ctq becomes higher again [Fig. 1], and 
therefore, Qhii quickly approaches to 1. 

Each curves of Fig. 6 show that there are two solutions to the equation Qhii{z) — 1, 
the first is at 2;^ > 15, and the second is at 2;^^ ~ 6 — 8. This corresponds to the two 
phases of the star formation. During the period between the redshift Zr and Zgp, the IGM 
is in multi-phase. The IGM consists of neutral and ionized hydrogen, HI and HII. Both 
components are comparable, and opaque for Lya photons. However, the temperature of 
such multiphase IGM can be of the order of 10^ K in average. Therefore, it is self-consistent 
to apply the bottom branch of the ctq of Fig. 1. 

When z ^ Zgp, the fiUing factor of neutral hydrogen becomes small quickly. The IGM 
is then transparent for Lya photons when the volume filling factor of HI is less than 10"^, 
and thus the redshift Zgp corresponds to the onset of the Gunn-Peterson transparency. 

It should be pointed out that the two peak curve of Qhii actually depends on the 
sharply dropping of ctq at Zr (Fig. 1). As has been discussed in §3.1, more realistic model 
should replace the sudden increase of IGM temperature by a soft process. In this case, Qhii 
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probably never reached 1 at Zr, but at most shows a peak around Zr- Therefore, before Zgp, 
the IGM has never been ionized to the level of 1 — Qhii < 10~^. Thus, the redshift Zgp is 
the first time of the IGM becoming the Gunn- Peterson transparency. Table 3 gives Zgp for 
three sets of parameters. 

Figure 6 and Table 3 shows that the redshift interval always as large as 

~ 10 — 12, which is not sensitive to the adopted parameters. This is because long duration 
from Zr to Zgp is essentially due to the zigzag z-dependence of the (Tq, which is less-sensitive 
to cosmological parameters, and the details of start formation. The temperature of the 
IGM increases by a factor of about 10^ during the reionization, while the entropy parameter 
T/riff^ also increases by few orders. Since high temperature and high entropy prevent 
baryonic gas from collapsing into the gravity well of dark matter, the reionization is slowed 
down with the increase of temperature and entropy of the IGM. It behaves like a negative 
feedback mechanism: once the rate of star formation is higher, the UV background produced 
by the star formation is also higher, and then more baryonic gas will be heated to higher 
temperature and higher entropy state. This finally yields a suppression on the rate of star 
formation. Therefore, the long duration of reionization is mainly given by the evolution of 
the temperature and entropy of the baryonic gas during the reionization. 

3.4. Optical Depth to the CMB 

If the reionization lasts the long period of Zr-Zgp, the number density of free electrons 
is significant through the whole epoch Zr > z > Zgp. One can then expect that in this 
reionization history the electron scattering optical depth to the CMB would be much larger 
than models with Zr — Zgp. With the solution of Qhii, the electron scattering optical depth 
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to the CMB can be calculated by (e.g. Hui et al. 2003) 



Te = 0.0525 X Qfc/i / 



VHn{z){l + zfdz 



(16) 



Jq 



where Vhii{z) is the volume filling factor of HII. Eq.(16) considered only the electrons 
from ionized hydrogen. To consider the contribution of He, we had assumed the filling 
factor of Hell to be identical to HII, i.e. the optical depth is increased by a factor of 
1.08. The Hell Helll ionization occurred at redshift ~ 3.3 (Theuns et al. 2002), and 
therefore, the correction of Helll on Tg is no more than 0.001. 

The relation between the volume filling factor Vhii and Qhii is dependent on the 
statistical behavior of the formation of the ionized region. If the ionized hydrogen region 
increase, dQuiiif) in the duration t to t -\- dt is without spatial memory, the formation 
probability would be uniform in space regardless whether there are ionized region formed 
at time before t. Thus, we have equation dYmi/dt = (1 — VHii)dQHii/dt, and the solution 
is Vhii{z) — 1 — exp[—QHii{z)]. However, the formation of the ionized region probably is 
not without spatial memory. The probability of forming collapsed clouds at the position 
where already hosted collapsed clouds should be low. Therefore, the process of the ionized 
hydrogen region formation most likely is self-avoid, i.e. the process has less probability 
at the positions where the process has been visited(e.g. Cardy 1996). For a perfectly 
self-avoiding process, no position can host the formation process more than once (e.g. 
Cardy 1996). In this case, we have Vhii{z) — Qhii{z) till Qhii < 1- 

Table 3 lists the results of Zr, Zgp and Tg for the model WMAP2 with parameters 
{Nicm,p/p) specified by (6, 40), (8, 40) and (6, 80). The optical depth Te(l) is for the model 
without any self-avoiding, while Te(2) is for model with perfectly self- avoiding. Real number 
probably is in between of these two ultimate models. All results shown in Tabic 3 give a 
low ~ 6 — 8 and a high Te~ 0.12 — 0.19. These results are excellently consistent with 
the observations of the redshift of the Gunn-Peterson trough and the electron optical depth 
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Table 3: Reionization redshift, Gunn-Peterson redshift and optical depth to CMB 



pIp 


N- 

^ ''ion 




Zgp 


re(l) 


re(2) 


6 


40 


18.0 


7.0 


0.150 


0.181 


8 


40 


16.0 


6.0 


0.123 


0.154 


6 


80 


20.4 


7.9 


0.159 


0.189 



to the CMB. The key of reconcihng the two observation is that there exists a long period 
of the universe to have 0.1 < Vhii < 1 — 10~^. The free electrons in this period contribute 
to a high value of Tg. On the other hand, the Gunn-Peterson transparency occurred at the 
end of this epoch, and the redshift Zgp is small. 



4. Discussion and Conclusions 

We showed with the LN model that the reionization and Gunn-Peterson transparency 
emerged at very different time. The thermal history of the universe around reionization can 
roughly be divided into three epochs: 1.) cold dark age z > Zri'va. which baryonic gas is cold, 
and opaque for Lya photons; 2.) hot dark age Zr > z > Zgp, in which a substantial part of 
baryon gas is ionized and hot, but it is still opaque for Lya photons, 3.) bright age z < Zgp. 
For the LCDM models favored by the WMAP, we have ~ 13 — 19 and 2;^^ ~ 6 — 8. That 
is, the hot dark age is in the redshift interval from 13-19 to 6 -8. This thermal history 
naturally yields a high optical depth = 0.12 — 0.19 and a low = 6 — 8, which is in 
excellent agreement with the observations of Tg (WMAP) and Zgp (Gunn-Peterson trough). 

The long duration from Zr to Zgp is essentially due to the evolution of the temperature 
and entropy of the baryonic gas during the reionization. For all the considered cosmological 
models, there is always a drop-off of ctq at Zr- Therefore, a large difference Zr — Zgp is a 
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common property of these models. The result of low Zgp and high Tg is less sensitive to the 
cosmological parameters. 

This result is also weakly dependent on the details of production of ionization photons 
from stars. This is due to the feedback mechanism of the start formation (§3.3). The higher 
rate of ionization photon production will cause a stronger input of entropy to baryonic gas. 
It then leads to a lower rate of the collapsing of baryonic clouds, suppressing the rate of 
star formation and the UV photon production. Therefore, model, which properly considers 
the dramatic evolution of the temperature and entropy of the IGM, will generally yield a 
low Zgp and a high Tg. 

The LN model provides an intuitive and quantitative description of the effect of the 
IGM temperature evolution upon the reionization history. This is because the LN model use 
directly variance ctq as a key parameter. It has already been recognized that the variance 
(Jo is more effective and successful in modehng the IGM detected by QSO's Lya forests 
(Bi & Davidsen 1997). We showed in this paper that the LN model provides a plausible 
description of photoionized and neutral IGM not only in the redshift range z < A, but 
extended to z ~ 15. 

The hot dark age between the redshifts Zr and Zgp is an interesting issue. During 
that period, star formation is significant, although the rate of star formation is lower than 
that at Zr- On the other hand, baryonic gas is in multiphase, and consists of comparable 
components of HI and HII. The gas has been already hot. Therefore, it is expected to find 
some observable signatures, such as high redshift 21 cm emission, of the hot dark age. 

JRL thanks Mr. Guopeng Shen and Jun Yu for useful discussions. LLF acknowledges 
support from the National Science Foundation of China (NSFC) and National Key Basic 
Research Science Foundation. 
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Fig. 1. — The variances of the hnear density fluctuations at the Jeans length as a function 
of redshift for model WMAP2. The reionization redshift z^. is taken to be 17. The dot line 
&X z < Zr is for the variances without considering the reionization. The solid line aX, z < Zr 
is considered the dropping of ctq at Zj. due to the increase of the temperature of IGM. 
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Fig. 2. — The cumulative mass fraction M(> p/ p) for model WMAP2. p/ p is taken to be 
2, 6 and 10. At 2; < Zr, the gq is taken from the dot line in Fig. 1. 




Fig. 3. — The clumping factor C(z) for the three cosmological models given in Table 1. 
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Fig. 4. — The filling factor Qhii vs. redshift of the three cosmological models, WMAPl, 
WMAP2 and COBE. p/p and Nion are taken to be 6 and 40, respectively. 




Fig. 5.— The filling factor Qhii vs. redshift of the model WMAP2, p/p = 6 and Ni, 
The bias parameter h of the clumping factor Chii — bC is taken to be 1.5, 2 and 4. 
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Fig. 6. — The redshift evolution of the volume fiUing factor QhiiCZj) from z > Zr to Zgp 
for model WMAP2. p/p and Nion are taken to be (8,40) (dot), (6,40) (soUd), and (6,80) 
(dash). The reionization redshifts Zr is given by Qhii ~ 1 at 2; > 10. The redshift of the 
Gunn-Peterson transparency Zgp is given by Qhii ~ 1 at z < 10. 



